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Itinerant electrons in a two-dimensional kagome lattice form a Dirac semi-metal, similar to 
graphene. When lattice and spin symmetries are broken by various periodic perturbations this 
semi-metal is shown to spawn interesting non-magnetic insulating phases. These include a two- 
dimensional topological insulator with a non-trivial Z2 invariant and robust gapless edge states, as 
well as dimerized and trimerized 'Kekule' insulators. The latter two are topologically trivial but the 
Kekule phase possesses a complex order parameter with fractionally charged vortex excitations. A 
charge density wave is shown to couple to the Dirac fermions as an effective axial gauge field. 



Certain physical observables in solids, such as the mag- 
netic flux in a superconductor or the Hall conductance in 
a quantum Hall liquid, are precisely quantized despite 
the fact that the host material may contain a significant 
amount of disorder. In all known cases this quantization 
phenomenon can be attributed to the notion of topologi- 
cal order. The bulk of such systems are characterized by a 
topological invariant that is insensitive to microscopic de- 
tails and robust with respect to weak disorder. Recently, 
a new class of topological invariants has been established 
to characterize all time-reversal (T) invariant band insu- 
lators in 2 and 3 spatial dimensions [U [21 [3l H] . These 
new invariants are of the Z 2 variety and the precisely 
quantized physical observable is the number of gapless 
edge (surface) states modulo 2. Topological insulators 
(TI) exhibit an odd number of edge (surface) states while 
trivial insulators exhibit an even number, possibly zero. 
In many ways the edge (surface) states of a TI behave 
as a perfect metal and are predicted to exhibit various 
unusual properties [6l |7] . They also show promise as 
possible components of future quantum computers [5] . 

Experimentally, HgTc/(Hg,Cd)Te quantum wells of 
certain width and composition have been identified as 2D 
topological 'spin Hall' insulators with robust edge states 
[§]. In addition, several 3D compounds involving bismuth 
have been so identified [101 E] an d several more have 
been predicted [121 H3] as likely candidates. In view of 
these rapid developments it appears likely that TIs might 
be a fairly common occurrence in nature. Given their ex- 
otic properties and their potential for technological ap- 
plications it is important to identify and study various 
model systems that exhibit this behavior. Such theoreti- 
cal understanding will aid experimental searches for new 
materials and help understand their unusual properties. 

In this Brief Report we advance the above agenda by 
describing a new class of 2-dimcnsional topological in- 
sulators on the kagome lattice, Fig. [TJu Although the 
properties of spin systems on the kagome lattice have 
been extensively studied, relatively little attention has 
been paid to the non-magnetic insulating phases of itin- 
erant electrons. In what follows we demonstrate, both 
analytically and numerically, that a simple tight-binding 
model of electrons on the kagome lattice at both | and 
I filling becomes a topological insulator upon inclusion 




FIG. 1: (Color online) a) kagome lattice is a triangular Bra- 
vais lattice with a 3-point basis forming a shaded triangle, b) 
The first Brillouin zone with the nodal points K± and time- 
reversal invariant momenta T n marked, c) The tight-binding 
band structure Eq. j3j. d) In the dimerized phase hopping 
amplitude along the thick (thin) bonds in ai direction is t + r/i 
(t—r)i). In the spin-orbit phase spin-up electrons hop between 
second neighbor sites with amplitude iX when moving along 
the arrow, —i\ against the arrow. For spin-down electrons 
the arrows are reversed, e) Trimerized 'Kekule' phase. 

of the spin-orbit (SO) coupling. Other T-invariant in- 
sulating phases include the dimerized and a trimerized 
band insulators; both have trivial Z2 invariants but the 
latter possesses a complex order parameter with vortices 
that carry fractional charge, in analogy to what happens 
in the Kekule phase of graphene [2]. We demonstrate 
that a charge density wave (CDW) modulation of the 
on-site energies does not produce a spectral gap but in- 
stead generates terms that couple as a U(l) gauge field 
to the low-energy Dirac fermions at | filling. Consistent 
with T-invariance, the gauge field couples with opposite 
sign to the two species of Dirac fermions, thus furnishing 
a concrete realization of an axial gauge field in a solid 
state system. 

We now proceed to substantiate the above claims. Our 
starting point is the tight-binding model 

H = -tJ2 C U- (1) 

where c\ a creates an electron with spin a on the site of 
the kagome lattice and (ij) denotes nearest neighbors. In 
momentum space Eq. (Il| becomes H = ^ kcr ^^'H^ka 
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with <J kCT = (ci kcr .c 2 k C r,C3k t T) T and 



T-/0 — 



cosfci COSfc 2 N 
— 2f I cosfci cosfc 3 

COS fc 2 COS k 3 



(2) 



The index Z = 1, 2, 3 in Qk<r labels the three basis sites in 
the triangular unit cell; ai = x, a 2 = (x + y3y) /2 and 
3-3 = &2 ^i denote the three nearest neighbor vectors, 
and — k ■ a^ . 

The spectrum of 7i£, Fig. lip, consists of one flat band 

(3) 

E^ — 2t and two dispersive bands 



E, 



(1,2) 



-1 ± y/AA* - 3 



(3) 



with Ak = cos 2 ki + cos 2 fc 2 + cos 2 £3. Bands 1 and 2 
touch at two inequivalent Dirac points K± = (±27r/3,0) 
located at the corners of the hexagonal Brillouin zone 
(BZ), Fig. [Tp>. Bands 2 and 3 touch at the center of the 
BZ. 

At I filing the lowest band is filled and the low-energy 
electronic excitations of H resemble those of graphcnc. 
To produce an insulator we now seek terms bi-linear in 
the electron operators that lead to the formation of a gap 
at the Dirac points. First, we focus on perturbations that 
do not further break the translational symmetry of H 
and preserve T. We have been able to identify two such 
terms: (i) a spin-independent lattice dimerization and 
(ii) a spin-orbit interaction induced hopping between the 
second neighbors [TJ H5J HE] • The former term breaks the 
inversion symmetry of the lattice while the latter breaks 
the SU(2) spin symmetry. 

The dimerization is described by -ffdim = 
^*ii c io- c J>' where 5t l0 = ±r] n describes an 
alternating pattern of bond hopping integrals along the 
three principal spatial directions as illustrated in Fig. 
[lji. In fc-space this becomes 



'0 



1_/dim 

tlx. 



irji sin k\ 




irjy sin fc 2 
1773 sin k 3 




(4) 



for both spin projections. (The lower triangle of the ma- 
trix is understood to be filled so that the matrix is hermi- 
tian.) The full expression for the spectrum of H.^ + 7i£ lm 
is complicated but it is easy to see that a gap Adi m — 2|ry| 
with 77 = 7/1 + rj2 + 773 exists at the Dirac points. 
The spin-orbit term takes form 



2\ 



d'ij) ■ (TapclaCjp, (5) 



where A is the spin-orbit coupling strength, d,-j are near- 
est neighbor vectors traversed between second neighbors 
i and j, and er is the vector of Pauli spin matrices. Since 

1 2 

djj all lie in the xy plane in our 2D model, only 0-3 
appears in Eq. ^ and the Hamiltonian decouples for 



the two spin projections along the z axis. The pattern 
of spin-orbit induced second neighbor hoppings then re- 
sembles the Haldane model [T3] and is illustrated in Fig. 
[T]i. In fc-space one obtains 



Hl° = ±2A 



icos(/c 2 




k 3 ) 



—icos (fc 3 
i cos (fci -4 




k 2 ) 
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where the +(— ) sign refers to spin up (down) electrons. 
Once again, although the full spectrum is complicated it 
is easy to deduce that a gap Aso = 4v3| A| opens up at 
the Dirac points. We remark that i?so opens a gap also 
between bands 2 and 3 while -f/dim does not. 

In order to develop some intuition for these insulating 
phases it is useful to examine the form of the low-energy 
Hamiltonians governing the excitations in the vicinity of 
the two Dirac points. This is obtained by linearizing 
H k = H^ + H^ im + Hi near K± and subsequently pro- 
jecting onto the subspace associated with bands 1 and 2. 
Assuming that both 7/ and A are non-zero we obtain 4 
independent Dirac Hamiltonians, 



^k£o 



v(r 3 k x + nky) + T 2 m ia , 



(7) 



labeled by the spin index a = ± and 'valley' index I = ±. 
Tj are Pauli matrices acting in the space spanned by the 
degenerate eigenstates of /V^ + and , v = v3t is the 
Fermi velocity and 



mi 



a = 2V3a\ + ir] 



(8) 



are Dirac masses whose relative signs define two distinct 
phases of the system. When \r)\ > 2y/3\X\ dimerization 
dominates and the Dirac masses at K + and K_ exhibit 
opposite signs, independent of spin. When \t}\ < 2v3|A| 
the SO interaction dominates and the mass signs at the 
two Dirac points are the same for a given spin but differ 
for the opposite spin projections. The two phases meet 
at a pair of critical lines | ^7 1 = 2\/3|A|. When crossing 
these lines, two out of four gaps close and the associated 
Dirac masses change signs. 

To understand the significance of the mass signs con- 
sider a boundary between the two phases, running along, 
say, the x = line in real space. For concreteness and 
simplicity we take A > 0, r\ — in the left half- plane and 
A = 0, rj > in the right half-plane. Focusing first on 
the Dirac point K + , we note that the spin-up mass m ++ 
remains positive for all x, suggesting a fully gapped spec- 
trum everywhere. The spin-down mass, however, neces- 
sarily undergoes a sign change across the x = boundary. 
Such a soliton mass profile is known to produce massless 
states [T7] in the associated Dirac equation, localized near 
the boundary. Specifically, Dirac equation 

[v(-ir 3 d x - irxdy) + T 2 m(x)]<p(x, y) = E<f>(x, y) (9) 

with m(x — > ±00) = ±toq has a gapless solution 



4>k(x,y) 



e lkv e 



US m{x')dx' 



(10) 
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extended along the boundary but localized in the trans- 
verse direction, with linearly dispersing energy Ek = vk. 
Similar analysis leads to a gapless state at node K_, but 
now for spin up and with Ek — —vk. 

A pair of spin- filtered, oppositely dispersing gapless 
edge states is a hallmark of the topological spin-Hall in- 
sulator [1, 2. 3 . The fact, apparent from the above con- 
struction, that these edge states depend only on the bulk 
band structure of the two insulators and not on the de- 
tails of the edge indicates their topological origin. It 
follows that one of these phases must be a topological 
insulator. It is easy to see that the dimerized phase is 
smoothly connected to a trivial insulator. Consider in- 
creasing all rj's continuously until r\\ = r/2 = r/3 = t. At 
this point the kagome lattice breaks down into a collec- 
tion of disconnected elementary triangles. This, clearly, 
is a trivial insulator. The spin-orbit phase, on the other 
hand, cannot be smoothly deformed into a trivial insula- 
tor and we show below by an explicit calculation that it 
indeed possesses a nontrivial Z2 invariant [TJ[21|S]. 

When a crystal possesses inversion symmetry the Z2 
topological invariant v is easy to evaluate. According to 
Ref. [TH] v is related to the parity eigenvalues £2™^) of 
the 2m-th occupied energy band at the four T-invariant 
momenta IV Our system is inversion symmetric when 
all 77's vanish and so we can use this method to find v. If 
we select site 1 of the unit cell as the center of inversion 
then the parity operator acts as V[tpi(r), ^(r), V'3( r )] = 
[ipi(— r), ip2{— r — 2ai), ip3(— r — 2a2)] on the triad of the 
electron wavefunctions in the unit cell labeled by vector 
r. In momentum space the parity operator becomes a di- 
agonal 3x3 matrix V k = diag(l, e" 2iai k , e - 2ja2 ' k ). The 
four T-invariant momenta in our system are marked in 
Fig. [ljj and can be expressed as Ti — n(x + y/V3)ni/2 + 
7r(— x + y/\/3)mi/2 with rij,mj = 0, 1. It is straightfor- 
ward to obtain the eigenstates of TLvi numerically and 
determine the parity eigenvalues of the occupied bands. 
We find that three £'s are positive and one is negative. 
Which of the four £'s is negative depends on the choice of 
the inversion center but the product ^£(1^) = (— l)" is 
independent of this choice and determines the non-trivial 
Z2 invariant v = 1, confirming our hypothesis that the 
spin-orbit phase at | filling is a topological insulator. 
Similar considerations for | filling also yield v = 1 . 

When the dimerization is present the inversion sym- 
metry is broken and we must use the more general 
method of Ref. [T] to find v. This relies on counting 
the number of pairs of first order zeros of the quantity 
P(k) = Pf (u m k|0|M„k), where « m k is the m-th eigen- 
state of 7ik, B is the time-reversal operator and the 
Pfaffian is taken over occupied bands m, n. We obtain 
timk by the numerical diagonalization of Tik and use it 
to straightforwardly evaluate P(k) . We find that the lat- 
ter contains exactly one pair of first order zeros (located 
at K ± ) in the SO phase when |?7| < 2\/3|A|, indicating 
v = 1. When |?7| > 2\/3|A| the zeros disappear and the 
system becomes a trivial insulator with v = 0. 

To further support our identification of the SO phase 




FIG. 2: Edge states in the lattice model for (a) spin-orbit 
and (b) dimerized insulator. A strip of width N y — 16 unit 
cells with open boundary conditions along y and infinite along 
x is used with \/t — 0.1 and ni/t — 0.1 for (a) and (b), 
respectively. 

as a topological insulator we have performed numerical 
diagonalizations of the lattice Hamiltonians Hq + i/dim 
and Hq + i/so using the strip geometry. In accord with 
the above arguments we find a pair of robust spin-filtered 
gapless states associated with each edge in the SO phase, 
Fig.|2ji, traversing the gap between bands 1 and 2. A sim- 
ilar pair of states traverses the gap between bands 2 and 
3, confirming that the SO phase is a topological insulator 
at I and | filling. In the dimerized phase there are gener- 
ically no such robust edge states, Fig. [SJd, although for 
certain values of parameters and types of edges gapless 
edge states can occur. The latter mirror the edge states 
found along certain types of edges in graphene [H] but 
are not topological in character. 

In graphene, a staggered on-site potential is known to 
open up a gap at the Dirac points [20]. We thus investi- 
gate the effect of the analogous CDW term, 

W£DW = diag(Mi;M2)M3)) (11) 

on the kagome semi-metal, where \xi represent the on- 
site potentials of the I = 1,2,3 basis sites of the lattice, 
independent of spin. Repeating the procedure leading to 
Eq. we find the low-energy Dirac Hamiltonians 

h kea = v [r 3 (k x - A e x ) + n(ky - A l y )] + tn, (12) 

for node £ and spin a, with (i — (/ii + /i2 + A*3)/3 and 

Ai = (2^3 - Mi - M2^/6\/3*, 

A e y = fa - iM,)i/6t. (13) 

The CDW couples to the Dirac fermions as a gauge field. 
Owing to the factor i in Eq. ( [T3| , the gauge potential has 
the opposite sign at the two Dirac points. Consequently, 
A must be thought of as an axial gauge field. This is 
consistent with the fact that the CDW does not break 
T. 

If we allow for additional fermion bi-linears that break 
the translational symmetry of Hq then many other in- 
sulating phases become possible. Of these, we briefly 
mention but one that realizes the analog of the Kekule 
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phase in graphene [2]. This occurs for a perturbation 
with a wave vector spanning the two nodal points K + and 
K_ in Fig. [TJ>. In the linearized theory such a pertur- 
bation leads to a Dirac mass that is off-diagonal in the 
space of nodal Hamiltonians and is in general complex- 
valued. In graphene, vortices in such a complex mass are 
known to carry a fractional charge ±e/2 (assuming spin- 
less electrons) [TJ] and obey fractional exchange statis- 
tics [21]. On the kagome lattice one possible realization 
of such a perturbation is the 'trimerization' depicted in 
Fig. [T£. The unit cell contains 9 atoms and there are 3 
distinct degenerate ground states related to the one de- 
picted by translations through 2ai and 2&2- We have 
verified that this perturbation indeed opens up a gap 
at the Dirac points and the three patterns produce off- 
diagonal Dirac masses with different complex phases. In 
analogy with the domain wall in the trimerized ID chain 
[22j we expect the vortex in the complex mass to bind 
fractional charge ±e/3, ±2e/3. This is also suggested 
by the Goldstone-Wilczek type counting argument [53] 
generalized to 2 spatial dimensions [24] that can be con- 
structed for the kagome lattice [55] . We leave the details 
of this and other interesting related problems, such as 
the exchange statistics of these objects, to future study. 

A question naturally arises as to the experimental real- 
ization of a system of itinerant electrons on the kagome 
lattice that would support the exotic phenomena pre- 
dicted above. Copper atoms in ZnCu3(OH)6Cl2, known 
as Herbertsmithite [26], and iron atoms in Jarosites [27] . 
provide good model Heisenberg kagome antiferromag- 



nets. The kagome lattice has also been argued to play 
a role in transport and magnetic properties of layered 
cobalt oxides, such as Na :c Co0 2 [28]. We hope that 
our theoretical findings will provide motivation for fu- 
ture experimental searches into possible realizations of 
the kagome lattice with itinerant electrons close to | or 
| filling. We note that it might be possible to artifi- 
cially engineer Hamiltonian by modulating the two- 
dimensional electron gas with a periodic potential with 
kagome symmetry, as recently demonstrated for 'artificial 
graphene' [2"9"] . 

When the basic tight-binding Hamiltonian is at hand 
the additional terms required to form an insulator can 
come about in various ways. In a crystal the SO cou- 
pling arises naturally. In graphene the relevant coupling 
strength is too small to open up a significant gap but in a 
lattice made of heavier ions A will be larger. The pertur- 
bations considered above can also arise as interaction- 
driven instabilities of itinerant electrons on the lattice. 
For the Dirac semi-metal at | filling a finite interaction 
strength U (of order t) is needed to open up a gap [30] . 
At | filling, however, the band crossing is quadratic and 
the instability towards the insulating phase occurs at in- 
finitesimal (repulsive) interaction |31j . making such sys- 
tem a promising candidate for a 2D topological insulator. 
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